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The precision limit in quantum state tomography is of great interest not only to practical applications 
but also to foundational studies. However, little is known about this subject in the multiparameter 
setting even theoretically due to the subtle information tradeoff among incompatible observables. 

In the case of a qubit, the theoretic precision limit was determined by Hayashi as well as Gill and 
Massar, but attaining the precision limit in experiments has remained a challenging task. Here we 
report the first experiment which achieves this precision limit in adaptive quantum state tomography 
on optical polarization qubits. The two-step adaptive strategy employed in our experiment is very 
easy to implement in practice. Yet it is surprisingly powerful in optimizing most figures of merit of 
practical interest. Our study may have significant implications for multiparameter quantum estima¬ 
tion problems, such as quantum metrology. Meanwhile, it may promote our understanding about 
the complementarity principle and uncertainty relations from the information theoretic perspective. 


I. INTRODUCTION 

Quantum state tomography is a procedure for infer¬ 
ring the state of a quantum system from quantum mea¬ 
surements and data processing [I|-[j|. It is a primitive 
of various quantum information processing tasks, such 
as quantum computation, communication, cryptography, 
and metrology 0-@|. In sharp contrast with the classi¬ 
cal world, any measurement on a generic quantum sys¬ 
tem necessarily induces a disturbance, limiting further 
attempts to extract information from the system. There¬ 
fore, many identically prepared systems are usually re¬ 
quired for reliable state determination. Conversely, the 
precision limit in quantum state tomography offers a per¬ 
fect window for understanding the distinction between 
quantum physics and classical physics IM1- 

Recently, great efforts have been directed to improv¬ 
ing the tomographic efficiency given limited quantum re¬ 
sources 1 13, 14], For example, adaptive measurements 
have been realized in experiments, which may improve 
the scaling of the infidelity in certain scenarios |14l . Il5l | . 
However, most studies have been tailored to deal with 
specific figures of merit under special settings, such as 
pure state or single parameter models, which admit no 
easy generalization to the more challenging and exciting 
multiparameter estimation problems with general figures 
of merit. In particular, the tomographic precision limit 
in the multiparameter setting is still poorly understood; 
experimental studies are especially rare. To fill this gap, 
in this work we report the first experiment that achieves 
the quantum precision limit in adaptive quantum state 
tomography on optical polarization qubits. 
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II. RESULTS 
Quantum Precision limit 

In practice, the tomographic precision limit is determined 
by experimental settings. As technology advances, it is 
ultimately limited by basic principles of quantum me¬ 
chanics. One fundamental limit is known as the quantum 
Cramer-Rao (CR) bound (UHlsj; see the appendix. In 
the one-parameter setting, this bound can be saturated 
locally by measuring a suitable observable, which may 
depend on the parameter point. To saturate the bound 
globally, it is usually necessary to employ adaptive mea¬ 
surements. A simple and effective choice is known as the 
two-step adaptive strategy |20L [2l| , whose basic idea can 
be sketched as follows. Suppose N copies of the true 
state are available for tomography. First, we perform 
a generic informationally complete measurement on N\ 
copies (usually Ni -C N especially when N is large) of the 
true state and compute the maximum likelihood estima¬ 
tor (MLE) |lj according to the measurement statistics. 
Then we perform the optimal measurement with respect 
to the estimator on the remaining = N — N\ copies 
and compute the MLE again. 

In the multi-parameter setting, however, the quantum 
CR bound generally cannot be saturated except when the 
optimal observables corresponding to different parame¬ 
ters can be chosen to be compatible. The existence of 
incompatible observables underlies the main distinction 
between quantum state estimation and classical state es¬ 
timation and is the main reason why multiparameter es¬ 
timation problems are difficult and poorly understood. 
Up to now, the optimal solutions have been found only 
for a few special cases 12, 1181, of which the qubit model 
is the most prominent |9j, l2ll . 1221 ] . 

To devise a good measurement scheme in the multipa¬ 
rameter setting, it is indispensable to take into account 
the subtle information trade-off among incompatible ob- 
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servables. A vivid manifestation of such tradeoff is the 
wave-particle duality relation between fringe visibility V 
and path distinguishability T> HH, 

V 2 + V 2 <1. (1) 

This phenomenon is not limited to the double-slit ex¬ 
periment but presents itself whenever we are trying to 
extract information about incompatible observables. It 
is especially important in understanding multiparameter 
estimation problems. Suppose the state of the quantum 
system is parametrized by a set of parameters denoted 
collectively by 6 , then such trade-off can be succinctly 
summarized bv the following inequality derived by Gill 
and Massar l2lil. 


tr{J- 1 (9)1 (9)} <d-l, (2) 

which is applicable to any measurement on a d-level sys¬ 
tem. Here 1(8) and J(9) are the Fisher and quantum 
Fisher information matrices, respectively (see the ap¬ 
pendix). The Gill-Massar (GM) inequality may be seen 
as a generalization of the wave-particle duality relation 
in the language of quantum estimation theory. To appre¬ 
ciate its significance, it is instructive to point out that 
the upper bound would be d 2 — 1 if all observables in 
quantum theory were compatible or, equivalently, if the 
quantum CR bound could always be saturated. 

The GM inequality imposes a fundamental precision 
limit on quantum state tomography based on individual 
(non-collective) measurements. For example, it sets a 
lower bound for the scaled weighted mean square error 
(WMSE) of any unbiased estimator [til. l2lj. 
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Figure 1. (color online) Quantum state tomography with two- 
step adaptive strategy. The observables cr' x , a' y , erf depend on 
the estimator pi obtained in the first step and are related to 
a x , a y , tj z via the unitary transformation U(s i) described in 
the main text. The probabilities pi, p 2 , and P 3 depend on 
both the estimator pi and the figure of merit. In the large-TV 
limit, it suffices to use the measurement statistics of step 2 
to construct the second MLE. In practice, it is preferable to 
employ the measurement statistics of both steps. 


as illustrated in Fig. [Q In each step, we need to imple¬ 
ment only three projective measurements that are mutu¬ 
ally unbiased. Despite the simplicity of this approach, it 
is capable of attaining the quantum precision limit with 
respect to most figures of merit of interest. To facilitate 
applications of our approach, we have determined the 
precision limits and local optimal measurements with re¬ 
spect to a large family of figures of merit in the appendix. 
For example, the GM bound for the scaled MSE is given 

by ana 


£ gm = (2 + Vl-s 2 ) 2 , (4) 

where s is the length of the Bloch vector of the true 
state. By contrast, the scaled MSE achievable by stan¬ 
dard tomography using mutually unbiased measurements 
is given by 3(3— s 2 ) [9|,[28|]. For concreteness, the two-step 
adaptive strategy for minimizing the MSE is sketched as 
follows: 


where W is the weighting matrix. Note that the GM 
bound for the WMSE is /N if the sample size is N. 
In the case of a qubit, the GM bound agrees with the 
bound derived by Hayashi [i22| and can always be satu¬ 
rated locally by mutually unbiased measurements [ij, [2l| 
(see the appendix). Recently, the GM inequality and 
GM bound have been turned into a powerful tool for 
studying a number of foundational issues entangled with 
incompatible observables 0 , such as the complementar¬ 
ity principle (25| and uncertainty relations [26|. The cross 
fertilization of quantum estimation theory and founda¬ 
tional studies is due to lead to deeper understanding of 
both subjects f?)l [~i~2l|. Determination of the quantum pre¬ 
cision limit in the multiparameter setting is thus of pri¬ 
mary interest from both practical and foundational per¬ 
spectives. 

Attaining quantum precision limit with two-step 
adaptive strategy 

Here we report the first experimental verification of the 
GM bound in adaptive quantum state tomography on 
optical polarization qubits. Our tomographic protocol 
consists of two-step adaptive measurements and MLE, 


1. Measure a x ,<j y ,a z on Ni/3 copies of the qubit, re¬ 
spectively, and compute the MLE p\ based on the 
measurement statistics. Denote the Bloch vector of 
h by s -i. 

2. Choose a unitary transformation U (si) that rotates 
Si to the z direction and apply this unitary trans¬ 
formation to the remaining N 2 = N — copies 
of the qubit state. Measure a x ,a y , a z with the fol¬ 
lowing probabilities (see the appendix for more de¬ 
tails), 


Pi =P 2 


1 

2 + y/1 - s 2 ’ 


P 3 


V 1 ~ Sj 
2 + y/1 -Sf’ 


( 5 ) 


where si is the length of si. Construct the MLE 
p 2 again based on the measurement statistics. 


Experimental setup 

The experimental setup is shown in Fig. [2] A 2-mm-long 
BBO crystal cut for type-I phase-matched spontaneous 
parametric down-conversion (SPDC) process is pumped 
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Figure 2. (color online) Experimental setup. A pair of hor¬ 
izontally polarized photons are generated via pumping a /3- 
barium borate (BBO) crystal. One is detected as a trigger 
and the other is sent through a half-wave plate (HWP), a 
quarter-wave plate (QWP), and a 770A quartz crystal in be¬ 
tween, marked as the state preparation module (green), to 
change the length and direction of the Bloch vector. After 
preparation, the photon polarization is measured by a con¬ 
figuration of QWP2, HWP2, and a polarizing beam splitter 
(PBS), marked as the adaptive measurement module (pink). 


at 404 nm by a 40-mW V-polarized beam. A pair of 
H-polarized photons with wave length A = 808 nm is cre¬ 
ated via SPDC process. One photon passes through an 
interference filter whose FWHM is 3 nm, resulting in a 
coherence length of 270A, and is detected by a single pho¬ 
ton detector acting as a trigger. The polarization state 
of the other photon is prepared by HWP1, a 770A quartz 
crystal, and QWP1. The quartz crystal, whose optic 
axis is aligned in the horizontal direction, completely de¬ 
coheres horizontal and vertical polarization components. 
The rotation angle of HWP1 determines the ratio of H 
and V polarization components, thereby together with 
the quartz crystal changing the purity of the state or the 
length of the Bloch vector. QWP1 is used to change the 
direction of the Bloch vector. QWP2 and HWP2 followed 
by a PBS are used to perform arbitrary projective mea¬ 
surements on the qubit. The rotation angles of HWP2 
and QWP2 are controlled by a labview program. The 
coincidence counts are measured by the coincidence cir¬ 
cuit and are analyzed by the computer. Without loss of 
generality, all polarization states in our experiments are 
prepared such that the Bloch vectors are aligned along 
the direction (0.490,-0.631,0.602). This is realized by 
setting QWP1 at the fixed deviation angle of 19.57°. The 
true state is calibrated with about 10 7 photons. 

To achieve high tomographic precision in experiments, 
it is crucial to reduce the systematic error to a very low 
level besides adopting the correct adaptive strategy. To 
this end, we have implemented error-compensation mea¬ 
surements proposed in Ref. (29| . wherein multiple nomi¬ 
nally equivalent measurement settings are applied to sub¬ 
ensembles such that main systematic errors cancel out in 
the first order. This method is capable of reducing the 
systematic error from 5 x 10 -5 to 8 x 10 -6 , which is 
about 100 times smaller than the statistical error when 
N = 9000 (see the appendix). 
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Figure 3. (color online) Precision limit with respect to the 
MSE. Experimental results of standard, adaptive, and known- 
state tomography are shown together with the theoretical 
MSE of the standard tomography and the GM bound. Here s 
is the length of the Bloch vector. In the experiment, N = 9000 
and Ni = 3000; each data point averages over 4000 repetitions 
and the error bar denotes the standard deviation of the mean. 
The MSE of the standard tomography is lower than the the¬ 
oretical value when the purity of the state is sufficiently high 
(depending on N) because the estimator is biased due to the 
influence of the boundary of the state space. 


Quantum precision limit with respect to various 
figures of merit 

To illustrate the power of the two-step adaptive strategy 
described above, here we investigate quantum precision 
limits with respect to a variety of figures of merit. In 
the first experiment, we verify this limit concerning the 
MSE with Ni = 3000 and N% = 6000. Here the values 
of N\ and N 2 are determined by numerical simulation to 
optimize the tomographic precision. The MSE is deter¬ 
mined by averaging over 4000 repetitions. As compari¬ 
son, we also implement two other tomographic strategies. 
In the first one, a Xl a y , a z are measured on iV/3 copies of 
the qubit, respectively, and the MLE is computed as be¬ 
fore. This standard tomography is widely used in prac¬ 
tice because of its simplicity. In the second one, a x ,a yi a z 
are measured with probabilities pi,P 2 ,P 3 as specified in 
Eq. ([5]) with si replaced by s after rotating the Bloch 
vector of each of the N copies of the qubit state to the z- 
axis. This “known-state tomography” assumes knowledge 
of the true state in designing the optimal measurement 
though not in reconstructing the state. It is not feasible 
in most practical situations, where the state is in fact 
unknown, but it is useful as a benchmark. 

Figure [3] shows the MSEs associated with the three 
measurement strategies mentioned above along with the 
theoretical MSE of the standard tomography and the 
GM bound. The experimental data agree very well with 
the theoretical prediction. In contrast with the adaptive 
strategy proposed in Ref. [1^, which offers no advan¬ 
tage over standard tomography with respect to the MSE, 
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Figure 4. (color online) Precision limit with respect to the 
MSB and WMSEs. In the experiment, N = 1200 and 
Ah = 300; each data point is an average over 1000 repeti¬ 
tions and the error bar denotes the standard deviation of the 
mean, (upper plot) MSBs of standard, adaptive, and known- 
state tomography together with the GM bound. The MSB 
of the known-state tomography is slightly smaller than the 
GM bound when the purity of the true state is sufficiently 
high (depending on N ) because the estimator is biased due to 
the influence of the boundary of the state space, (lower plot) 
WMSEs with respect to the family of monotone Riemannian 
metrics determined by Eq. © for a state with s = 0.9. 


our two-step adaptive strategy is significantly more effi¬ 
cient than standard tomography. What is remarkable, 
the MSE achievable by the two-step adaptive scheme is 
quite close to that of known-state tomography and satu¬ 
rates the GM bound approximately. 

Next, we investigate the quantum precision limit with 
respect to the mean square Bures distance (MSB). Inci¬ 
dentally, the GM bound for the MSB is equal to 9/(4 N) 
(see the appendix), which is independent of the unknown 
qubit state. The two-step adaptive strategy for mini¬ 
mizing the MSB is quite similar to that for minimizing 
the MSE except that in the second step <r Xl a y , a z are 
measured with probability 1/3 each. Similar modifica¬ 
tion applies to the “known-state tomography”. Unlike 
the previous case, it is not necessary to implement error- 
compensation measurements [29| because a simple mea¬ 


surement strategy can already achieve a sufficiently high 
precision. Here our two-step adaptive strategy is sim¬ 
ilar to the one proposed in Ref. [Tj], which is tailored 
to minimize the infidelity. This is not surprising since 
the infinitesimal square Bures distance is equal to the in¬ 
fidelity. Our study provides further justification of the 
approach in Ref. [14| from a wider context. 


The upper plot of Fig. 0] shows the MSBs associated 
with three measurement strategies along with the theo¬ 
retical precision limit set by the GM bound. In standard 
tomography, the MSB increases rapidly as the purity of 
the true state increases. In sharp contrast, the MSB 
achievable by the two-step adaptive scheme is almost in¬ 
dependent of the purity and is much smaller than that 
of standard tomography, the more so the higher purity 
of the true state. Moreover, it saturates the GM bound 
approximately, except when the true state is nearly pure. 
The small gap from the bound is mainly due to the 
fact that the experiment under investigation is not close 
enough to the asymptotic regime as the ratio Ni/N is 
nonnegligible. Meanwhile, the efficiency gap of the stan¬ 
dard measurement from the local optimal measurement 
is significant for states with high purities. 


A prominent merit of our approach is its versatility 
in dealing with various figures of merit as emphasized 
before. To further demonstrate this point, we now turn 
to verifying the precision limit with respect to WMSEs 
based on monotone Riemannian metrics [30l - l32j |: see the 
appendix. For concreteness, we shall focus on the family 
of metrics characterized by the following equation 


2 1 / ds 2 s 2 dfl 2 \ 

= lVT^ + ^(ifi)(IT7)r 


( 6 ) 


where dfl 2 is the metric on the unit sphere and / n (f) = 
[(1 + t 1 /")/ 2] n . The GM bound for the scaled WMSE 
turns out to be 


cGM 

Wn 


1 

4 


( 2 + iY h 

r(i + s) 1 / n + (i-s) 1 / n i 

\h + 1 ) ’ h ~ 

2 J 


( 7 ) 


The metric reduces to the Bures metric when n = 1 and 
the quantum Chernoff metric (33| when n = 2. Other 
metrics can also be treated on the same footing. The two- 
step adaptive strategy is quite similar to the one in the 
previous case, except that in the second step a x ,a Vl a z 
should be measured with the probabilities pi = P 2 = 
1/(2 + h) and P 3 = h /(2 + h), where h is defined in the 
same way as h in Eq. © but with s replaced by si. 
The lower plot of Fig. [4] shows the experimental results 
for a true state whose Bloch vector has length 0.9. The 
WMSEs achieved by the two-step adaptive strategy are 
much smaller than that of standard tomography as in the 
previous case. Moreover, they agree very well with GM 
bounds. 
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III. DISCUSSION 

We have implemented optimal adaptive qubit state to¬ 
mography of mixed states in the multiparameter setting. 
The two-step adaptive strategy employed in our exper¬ 
iment is very easy to realize in practice. Yet it is sur¬ 
prisingly powerful: It is applicable to optimizing most 
figures of merit of practical interest, such as the WMSE 
based on the Bures metric, quantum Chernoff metric, or 
any other metric commonly used. Moreover, it is capa¬ 
ble of attaining the precision limit set by the GM bound 
approximately. Our experiment represents a significant 
step towards optimal quantum state tomography in the 
multiparameter setting, which may have profound im¬ 
plications for multiparameter quantum metrology. Fur¬ 
thermore, our study exemplifies the subtle connection be¬ 
tween the complementarity principle and quantum pre¬ 
cision limit, thereby promoting the cross fertilization of 
quantum estimation theory and foundational studies. 
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IV. METHODS 


Data collection 

To stabilize the collecting efficiency, multimode fibers are 
used to direct photons from free space to single-photon 
detectors. The whole experimental setup is covered by 
a black hood to keep away stray light from outside; the 
multimode fibers are wrapped up by black plastic films to 
suppress laser scattering. To balance the data-collection 
time and random coincidences, the power of the contin¬ 
uous pumping laser is adjusted to 40 mW, and the coin¬ 
cidence window is set at 2 ns. The resulting coincidence 
rate is 8000 per second, and the coincidence efficiency is 
about 8%. 
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Appendix A: Quantum Cramer Rao bound 


Suppose the state p(9) of a given quantum system is 
characterized by a set of parameters 61 , 62 ,... , 8 g . To de¬ 
termine the values of these parameters, we may perform 
measurements and construct an estimator based on the 
measurement statistics. Once a measurement with out¬ 
comes IIj (II{ are positive operators that sum up to the 
identity) is chosen, the probability of obtaining each out¬ 
come is determined by the Born rule p(f\9) = tr(pllj). 
Given an estimator 0(£), the accuracy is usually quanti¬ 
fied by the mean-square-error (MSE) matrix, with entries 
given by 


C 3 k { 6 ) = E[(0j(£) - ^0(4(0 - e k )]. (Al) 


The famous Cramer-Rao (CR) bound in statistical infer¬ 
ence states that the MSE matrix C(9) of any unbiased 
estimator is lower bounded by the inverse of the Fisher 
information matrix I(9)\ that is, C(9) > I~ 1 (9), where 


Ijk( 6 ) = E 


d}npm\(d±pmx 
d6 :i )y de k ) ■ 


(A2) 


Meanwhile, I( 6 ) is upper bounded by the quantum Fisher 
information matrix J(9) with 


JAO) = \t*{miL j ( 8 )Lk{d)+L k (d)L j { 8 )\}, (A3) 


where Lj(9) is the Hermitian operator that satisfies the 
equation 


^. = l[p( 6 )L j ( 6 ) + L j ( 6 )p( 6 )} (A4) 

and is known as the symmetric logarithmic derivative 
(SLD) of p(9) with respect to 8j. Together with the clas¬ 
sical CR bound, the inequality 1(6) < J(6 ) implies the 
quantum CR bound C(9) > J~ 1 (9) prj - flaj . In the one- 
parameter setting, this bound can be saturated locally 
by measuring the observable L(9). In the multiparam¬ 
eter setting, however, it generally cannot be saturated 
except when the Lj(9) commute with each other. 


Appendix B: Gill Massar bound for the qubit 


Here we show that the CM bound for the WMSE in 
the qubit case can always be saturated, and the local 
optimal measurement can be realized by a complete set 
of mutually unbiased measurements S 0,0- Recall 
that the CM bound for the scaled WMSE corresponding 
to the weighting matrix W is given by 


cGM 

°W 


(tr v j i / 2 if j -1 / 2 y 


(Bl) 


This bound can be saturated if and only if there exists a 
measurement that yields the Fisher information matrix 


( B 2) 

When W and J commute, Eqs. m and (IB2I) reduce to 


gm _ (r vwr^y 

c w ~ 


d- 1 


(,- 1 yw 

tr VWJ^ 


In the case of a qubit, it is convenient to parametrize 
the state space by the components of the Bloch vector 
s = (s x , s y , s z ). Then the inverse quantum Fisher in¬ 
formation matrix takes on the form J _ 1 (s) = 1 — ss. 
Suppose that I\y in Eq. (IB2I) has eigenvalues 01 , 02,03 
along with orthonormal eigenvectors ri,r 2 ,r 3 . Denote 
by Si,S 2 ,S 3 the three components of the Bloch vector 
in this basis; denote by <x the vector composed of three 
Pauli matrices cr x ,a y ,a z . Then the GM bound can be 
saturated by measuring each observable ctj := rj ■ cr with 
probability a ? (l — s'f). The equality Oj(l — s 2 ) = 
tr (J^Iw) = 1 0 [llj guarantees that the set of proba¬ 
bilities is indeed normalized. This observation confirms 
our claim that the optimal measurement scheme can be 
realized using a complete set of mutually unbiased mea¬ 
surements. 

Next, we derive the GM bound when the weighting 
matrix is the identity or determined by a monotone Rie- 
mannian metric |30l I32l |. Without loss of generality, we 
may assume that the Bloch vector is aligned with the 
z-axis, so that J = diag(l, 1, 1/(1 — s 2 )), where s is the 
length of the Bloch vector. If the weighting matrix is di¬ 
agonal, say W = diag(uq, W 2 , W 3 ), then Eq. (IB3I) reduces 
to 


£w M = {VA + VA + VmV -s 2 )) 2 , 

diag(- v /uii, ,Jw 2 , y/w 3 (l - s 2 )" 1 ) ( B4 ) 

±w = - . —. 

y/wi + \A^2 + V W 3 (l - s 2 ) 

The optimal measurement scheme can be realized by 
measuring <jj for j = 1,2,3 (cr x ,a y ,a z in this case) with 
probability 


Pi = 


Vwjjl- s 2 S j3 ) 

y/m + ifwf + i/w 3 (l - s 2 ) 


(B5) 


The GM bound for the MSE is obtained when w\ = W 2 = 
IU 3 = 1, in which case we have (cf. Refs. S0) 


£ GM = (2+y/l -s 2 ) 2 , 


Pi = P2 


1 

2 + VTA^' 


P3 


VTA? ( B6 ) 
2 + VTA?' 


Every monotone Riemannian metric for a qubit has the 


d- 1 
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following form up to scaling, 

2 1 / ds 2 s 2 dfl 2 \ 

= 4 \T^ + f(±^)(l + s ))' 


(B7) 


where dfl 2 is the metric on the unit sphere, and / is 
a Morozova-Chentsov function [30l - [32| . If the weighting 
matrix is determined by this metric, that is, 


w 1 =w 2 = -[(1 + s)f(t)] \ 


U>3 = 


4(1 -s 2 )’ 


(B 8 ) 


from six aspects: uncertainties in the phases and rota¬ 
tion angles of the optical axes of QWP2 and HWP2, the 
extinction ratio of the PBS, and unbalanced collecting 
efficiencies between the two branches from the PBS. In 
the analysis of systematic errors, we neglect the statisti¬ 
cal fluctuation of measurement results. When the Pauli 
operator r ■ cr with unit vector r is measured on N copies 
of the qubit state p = (1 + s ■ cr)/2 , the expected photon 
counts of the outcomes ±1 are 

N± = Np±r/±, (Cl) 


where t = (1 — s)/(l + s), then 
1 

Pl=P2 = 


2 + \/(l + s)f(t) 


P 3 = 


Sf M = \{2((1+ s)f(t)}~ 1 / 2 + !} 2 . 


\/(l + s )f(t) 

2 + \/(l + s)f{t) 
(B9) 


An important family of Morozova-Chentsov functions has 
the form dH 


fn(t) = ( 


1 + 1 V" 


(BIO) 


Accordingly, Eq. (IB9I) reduces to 


1 


Pi = P2 = 


2 + h' 


P3 = 


2 + h’ 


where 


cGM _ 


h = 


1(2 A 

l(h + 1 ) 5 

(l + s) 1 /" + (l-s) 1/n W/2 


(Bll) 


(B12) 


where p± are probabilities given by the Born rule and ij± 
are collecting efficiencies determined by channel losses 
as well as coupling and detection efficiencies of the two 
branches. 


Consider a PBS with an extinction ratio of 1//3, the 
measurement operators corresponding to the two out¬ 
comes ±1 actually realized are II± = [1±(1 — 2/3)r -er\/2. 
Accordingly, p± = [1 ± (1 — 2/3)s • r\/2. The unbalance 
between collecting efficiencies of the two branches is de¬ 
fined as 77 = (t 7 +/? 7 _) — 1. To the first order in 77 , the 
frequencies of measurement outcomes are 

f± = n+ N +n- ^ p±± (C2) 

The expectation value of r ■ cr is 

to = /+ — /_« m — 2m/3 + 2p + p_r), (C3) 

where 

m = s ■ r (C4) 


The metric turns out to be the Bures metric when n = 1, 
in which case p\ = P 2 = P 3 = 1/3 and = 9/4, and 
the quantum Chernoff metric [33| when n = 2. 


Appendix C: Systematic error 


is the expectation value if the extinction ratio is infinity 
and the collecting efficiencies are balanced. In our exper¬ 
iment, the measurement of r ■ cr is realized by a config¬ 
uration of QWP2 and HWP2 as specified by the phases 
Si, 62 and rotation angles 9\, 82 of QWP2 and HWP2; see 
Table U for examples. 


In the ideal scenario, the error between the calibrated 
state and the true state vanishes in the large sample limit, 
so the calibrated state can serve as the true state. In real 
experiments, however, the error does not vanish due to 
various experimental imperfections. By the same token, 
different calibration procedures may result in different 
calibrated states. To realize the goal of our experiment, 
it is crucial that the total systematic error is much smaller 
than the error due to statistical fluctuation. This is espe¬ 
cially the case in the experiment concerning the MSE, in 
which the statistical fluctuation is suppressed by a large 
sample and large number of repetitions. In this appendix, 
we estimate the magnitude of the systematic error and 
show that it meets the requirement of our experiment. 

The systematic errors in our experiment mainly come 


According to Eqs. (IC3II and (IC4I) . the systematic error 
of to can be estimated as follows, 

(Ato ) 2 = £(to c ) 2 (AC) 2 , (C5) 

C 

where to q is the partial derivative of m with respect to £ 
with ( = /3,r], 61 , 62 , 81,82 taken at ideal values of these 
parameters, that is, /3 = 77 = 0 , 61 = ir/2, 62 = 7 r (the 
ideal values of 81 and 82 depend on r). These deriva¬ 
tives are (to^) 2 = 4to 2 , (to ,,) 2 = 4(p + p _) 2 < if 
C = <5i, 82 , 0i, 02, then (to c ) 2 = (me ) 2 = (s • r ) 2 < |r f | 2 , 
where = dr/d) Although the dependence of r on 
61 , 62 , 81,82 is complicated |29|, the magnitudes of its par- 

















Table I. Rotation angles of QWP2 and HWP2 for realizing 
the measurement of the Pauli operator r ■ a. Also listed are 
the partial derivatives of r with respect to phases and rotation 
angles of the two wave plates. 


r 

(. 01 , 62 ) 

|R5il 

\ r 5 2 \' 

kej 

\ r e 2 \ z 

(1, 0, 0) 

(45°, 22.5°) 

0 

0.5 

8 

16 

(0, 1, 0) 

(0°, 22.5°) 

1 

0.5 

4 

16 

(0, 0, 1) 

O 

O 

O 

0 

0 

8 

16 


tial derivatives have simple forms: 

|r 5l | 2 = sin 2 (20i - 40 2 ), 

K| 2 = 4 + 4cos 2 (20 1 - 40 2 ), 

\rg 2 \ 2 = 16. 

When r is one of the basis vectors, the magnitudes of 
these partial derivatives are listed in Table Q] 

In our calibration process, three Pauli operators r J ■ <x 
are measured, where r 0 for j = 1,2,3 form an orthonor¬ 
mal basis, so that the eigenbases of • cr are mutually 
unbiased. The calibrated state is determined by the for¬ 
mula 


P = 


1 + E,=i a 


(C7) 


The systematic error between the calibrated state and 
the true state can be calculated as follows, 


£ E(^c) 2 (AC) 2 = £(M c ) 2 ( AO 2 , (C8) 
C J=1 c 


where (M c ) 2 = Ej=i("4) 2 and C = /3,V,Si,d 2 ,0 3 1 ,9 3 2 . 

In order to decrease the systematic error caused by 
the PBS, we use a beam displacer (BD) acting as a 
PBS, whose extinction ratio is about 8000. The result¬ 
ing contribution to the systematic error is (Ms) 2 (AB) 2 < 
6 x 10 -8 . Notice that 


( M p ) 2 = £(^) 2 = 4£(rrP) 2 = 4|s| 2 < 4 (C9) 


1=1 


1=1 


since r J for j = 1,2,3 are orthonormal. 

The unbalance between the collecting efficiencies of the 
two branches is mainly due to the fluctuation caused by 
the disturbance of the air flow as well as the vibration 
and shift of the laser propagation direction. With a 
multimode fiber collecting system, this fluctuation can 
be reduced to below 0.002 within 16 hours, the mea¬ 


surement duration of 4000 repetitions for one quantum 
state. This fluctuation results in ( At f ) 2 = 4 x 10 -6 and 
(M r] ) 2 {Ari ) 2 < 3 x 10 -6 . Notice that (m^) 2 < \ and 

W = E- = 1 (K) 2 <!- 

The wave plates provided by our manufacturer have 
phase uncertainties of about |A<5| = 1.2° from ideal 
phases, that is, <5i = 90° ± 1.2° and S 2 = 180° ± 1.2° 
for QWP2 and HWP2. The phase uncertainties of QWP2 
and HWP2 together contribute [(Ms ^ 2 + (Ms 2 ) 2 ](AS ) 2 < 
8.8 x 10~ 4 to the systematic error. This systematic error 
would wash out the advantage of adaptive state tomogra¬ 
phy over standard tomography when N = 1200. To solve 
this problem, we use calibrated phases of <5i = 88.7°±0.3° 
and S 2 = 179.9° ± 0.3° provided by the manufacturer to 
calculate the rotation angles required to realize desired 
measurement settings. Consequently, the systematic er¬ 
ror is reduced to [(M^J 2 + (Ms 2 ) 2 ](AS ) 2 < 5.5 x 10 -5 . 

The rotation stages used in our experiments have a 
precision of 0.01°; the uncertainties in rotation angles 6 i 
and 02 of QWP2 and HWP2 are mainly determined by 
the uncertainties in the calibration angles of their optical 
axes, which are about A0 = 0.1°. The resulting system¬ 
atic error is [( Mg .J 2 + (Mg 2 ) 2 ](A9 ) 2 < 2 x 10 -4 . 

The total systematic errors due to the six aspects of 
experimental imperfections mentioned above is no larger 
than 2.6 x 10 -4 . Incidentally, the error between states 
calibrated by standard tomography and known-state to¬ 
mography with about 10 7 photons in the experiment is 
around 5 x 10~ 5 . In our experiments concerning the MSB 
and WMSEs with respect to a family monotone Rieman- 
nian metrics, this systematic is low enough to guarantee 
that the experimental result agrees well with the theoret¬ 
ical prediction and that the advantage of adaptive state 
tomography over standard tomography is clearly mani¬ 
fested. 

In the experiment concerning the MSE, however, the 
systematic error is still not small enough since the ad¬ 
vantage of adaptive strategy over standard strategy is 
quite limited and we need to adopt a large sample size 
and a large number of repetitions to overcome this is¬ 
sue. To solve this problem, we developed and imple¬ 
mented error-compensation measurements |2?)l | . wherein 
multiple nominally equivalent measurement settings are 
applied to sub-ensembles such that main systematic er¬ 
rors cancel out in the first order. This method alleviates 
the influences of uncertainties in the phases and the cal¬ 
ibration angles of QWP2 and HWP2. In this way, the 
total systematic error can be reduced to 7 x 10~ 6 . Ac¬ 
cordingly, the error between the two calibrated states by 
standard tomography and known-state tomography with 
about 10' photons is reduced to 8 x 10 -6 |2!|. Now the 
systematic error is much smaller than the MSE gap be¬ 
tween adaptive tomography and standard tomography. 
Consequently, the experimental result agrees well with 
the theoretical prediction again. 
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